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ABSTRACT 

We present a new time-dependent inhomogeneous jet model of non-thermal blazar 
emission, which reproduces the entire spectral energy distribution together with the 
rapid gamma-ray variability. Ultra-relativistic leptons arc injected at the base of a 
jet and propagate along the jet structure. We assume continuous reacceleration and 
cooling, producing a relativistic quasi-maxwellian (or "pile-up") particle energy distri- 
bution. The synchrotron and Synchrotron-Self Compton jet emissivity are computed 
at each altitude. Klein-Nishina effects as well as intrinsic gamma-gamma absorption 
are included in the computation. Due to the pair production optical depth, consider- 
able particle density enhancement can occur, particularly during flaring states. Time- 
dependent jet emission can be computed by varying the particle injection, but due to 
the sensitivity of pair production process, only small variations of the injected density 
are required during the flares. The stratification of the jet emission, together with a 
pile-up distribution, allows significantly lower bulk Lorentz factors, compared to one- 
zone models. Applying this model to the case of PKS 2155-304 and its big TcV flare 
observed in 2006, we can reproduce simultaneously the average broad band spectrum 
of this source as well as the TeV spectra and TeV light curve of the flare with bulk 
Lorentz factor lower than 15. 

Key words: galaxies: active - BL Lacertae objects: individual: PKS 2155-304 - 
galaxies: jets - gamma-rays: theory - radiation mechanisms: nonthermal 
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1 INTRODUCTION 

It is widely admitted that the blazar phenomenon is due to 
relativistic Doppler boosting of the non-thermal jet emission 
taking place in radio-loud Active Galactic Nuclei (AGN) 
whose jet axis is closely aligned with the observer's line of 
sight. Blazars exhibit very broad spectral energy distribu- 
tions (SED) ranging from the radio to the gamma-ray band 
and dominated by two broad band components. In the Syn- 
chrotron Self Compton scenario (SSC), the lowest energy 
hump is attributed to the synchrotron emission of relativis- 
tic leptonic particles, and the highest one is attributed to 
the Inverse Compton process (IC) of the same leptons on 
the synchrotron photon field. Broad band observations of 
these objects are crucial to understand the jet physics and 
to put reliable constraints on jet parameters. 

The most extreme class of blazars are the highly peaked 
BL lac sources (HBL), where the synchrotron/Inverse 
Compton components peak in the UV/X-ray /gamma-ray 
(GeV up to TeV) range. The recent development of obser- 
vational techniques in the TeV range, with Cherenkov tele- 
copes experiments like HESS, MAGIC or VERITAS has al- 
lowed the detection of about 18 HBL above 300 GeV. These 



objects are well known to be highly variable in all energy 
bands, from radio to gamma-ray, with timescales varying 
with energy. Perhaps the most extreme and remarkable ex- 
ample of this extraordinary variability behaviour has been 
caught by the HESS instrument in th e recent observations o f 
PKS 2155-304 during summer 2006 (|Aharonian et al.ll2007l ) 
(impressive observations has als o been recently obtained for 
Mkn 501 i.e. [Albert eralll2007l ). 

The most simple models of high energy emission assume 
a one-zone, homogeneous region. The SSC emission is as- 
sumed to come from a spherical zone of radius R, filled with 
relativistic leptons characterized by their density and a 
characteristic Lorentz factor 7 (e.g. that most contributing 
to the peak emission) , embedded in a isotropic magnetic field 
B. The blob is assume to move with a bulk Lorentz factor Tb 
yielding to a Doppler factor Sb — [^b{i — f3cos6)]~^ where 
6 is the angle between the blob direction of motion and the 
line of sight. In this model, all physical quantities are aver- 
aged on the whole spherical region, and, more importantly, 
the synchrotron and IC emission are cospatial. Causality 
constrains and a necessarily low pair creation optical depth 
then generally implies rela tively high bulk Lorentz factor 
iMastichiadis fc Kirk|[l997l) . 
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In the case of th e 2006 big flare of PKS 2155-304 
l|Aharonian et alll2007l ). the observed variabihty time scale 
(~ 200 sec) in the TeV range implies a minimum bulk 
Loren tz factor greater than 50 (jBegelman. Fabian, fc ReesI 
I2OO8I ) assuming an homogeneous one zone model. However, 
such high values of the bulk Lorentz factor are in contradic- 
tion w ith constrains derived f rom other observatio nal evi- 
dence (|Urrv fc Padovanilll995l : iHenri fc Saugil2006l and ref- 
erences therein). Furthermore, one-zone models are unable 
to flt the entire spectrum, the low energy radio points being 
generally attributed to more distant emitting regions. More 
complex models h ave been proposed including for example 
jet s tratification (|Ghisellini et al. 19851: Katarzviiski et al 



2005 ) or jet deceleration (e.g. iGeorganopoulos fc Kazanas 
20031 ). We present here a new approach, unifying small and 
large scales emission region: we consider that the radio jet is 
actually filled by the same particles originating from the high 
energy emitting region, at the bottom of the jet, that have 
propagated along it. We describe thus the emitting plasma 
by a continuous (although variable) particle injection, sub- 
mitted to continuous reacceleration and radiative cooling. 
This model fits well i n to th e two -flow framework originally 
proposed bvlPelletied (1 19851') and ISol et al.1 l\198^ (see also 
iTsinganos fc Bogovalovll200S or th e "spine-in-jet" model de- 
velopped bv lChiaberge et al.l2000l ) where a non relativistic, 
but powerful MHD jet launched by the accretion disk, sur- 
rounds a highly relativistic plasma of electron-positron pairs 
propagating along its axis. The MHD jet plays the role of 
a coUimater and an energy reservoir for the pair plasma, 
which is responsible for the observed broad band emission. 
The present model only concentrates on the physical pa- 
rameters of the relativistic pair beam and is described in 
Sect. [21 We show its application to the case of PKS 2155- 
304 in Sect. [3] focusing on the 2006 big flare event. We then 
discuss our results in Sect. |4] 



2 DESCRIPTION OF THE MODEL 
2.1 Geometry of the model 

We consider that the relativistic plasma propagates in a sta- 
tionary funnel whose geometry is parametrized as follows: 



r{z) = Ro 



i/c 



(1) 



where r(z) is the radius of the jet at the altitude 2. This 
shape describes a jet with a "shifted" paraboloid shape, with 
an initial inner radius Ri at z = 0, and a radius Ro at a 
distance Zq from the apex. The index w is lower than 1 for 
a coUimated jet. The magnetic held inside the jet is radially 
averaged at each altitude, and has the following dependence: 



(2) 



We take A between 1 and 2, the two extreme values describ- 
ing respectively a pure toroidal or a purely poloidal held 
distribution. 

We consider that the jet is continuously accelerating, 
starting from rest (Fj, — 1 nt z = 0) and reaching an asymp- 
totic value Fboo. The acceleration is assumed to take place 
over a distance comparable with Zo- A detailed model of the 
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Figure 1. Sketch of the jet geometry. See text for the signification 
of the diff'eront parameters. 



acceleration mechanism is beyond the scope of this paper, so 
we chose a simple parametrization to describe the evolution 
of Tbiz): 

I "1 l/a 



n{z) 



1 + 



ra 



1 + ^ 



(3) 



a being a "stiffness" parameter describing the width of the 
accelerating region. 



2.2 Particle energy distribution 

In our model, the energy distribution function (EDF) of 
the electron-positron plasma is assumed to be a relativis- 
tic maxwellian (or "pile- up") distribution : 



n{^,z,t) = no{z,t)'y exp 



7 



70(2, i) 



(4) 



where no{z,t) is a normalization factor and 70 the 
characteristic "pile-up" Lorentz factor. This distribu- 
tion is a natural outcome of some acceleration pro- 
cesses like sec ond or der Fermi acceleration or magnetic 
reconnection ([Henri 
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Sauge 



Pelletieij Il99ll : ISchlickeiseJ 



2004h 



. . - - - - ■ ■ It has been also shown by 

Henril ( 2004 ) that a pile up distribution is well 
suited to reproduce the narrow peaked high energy com- 
ponent of TeV blazars. 

When the plasma propagates in the structure, the parti- 
cles loose energy via synchrotron and Inverse Compton cool- 
ing, producing the observed emission. It turns out however 
that the cooling time is much too short for the relativistic 
particles to flU the whole jet. It is thus necessary to assume 
a fast reacceleration process along the jet. The acceleration 
rate is parametrized by a shifted power law with an index 
^ and an exponential cut-off after some altitude Zc to avoid 
energy divergence: 



Qacc{z) = Qo 



z 


( Ri \ 




\Ro) 



-c 



exp 



(5) 



In the two-flow framework (see introduction), this re-heating 
is naturally provided by the surrounding MHD structure via 
second order Fermi process. Then we only need to determine 
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Table 1. Model parameters of the flaring and quiescent state. During the flaring state, the flux of injected particles varies between the 
indicated minimum and maximum values, following the injection pattern displayed in Fig. 3. The other parameters remain fixed. Ri, Rq, 
Zq, Zc are in unit of 10^^ cm 



STATE ^(ZO NiZi) $(Zo) 7V(Zo) Qo R, Ro Zo Zc Bo cu X C 

[10*2s-l] [cm-3] [1042g-i] [cm-^] [s^i] [G] 

Max 2.09 1817 70.1 1300 

flaring Aver. 1.84 1666 24.4 600 

Min 1.16 1055 2^33 60 

quiescent TIE 1051 135 ST 



6.5 



2.5 



15 1.1 1.78 20 5 X 10' 



0.2 1.9 1.27 



two parameters to fully describe the relativistic plasma: 
no{z,t) and jo{z,t); 70(2;, t) is determined by balancing ra- 
diative losses and re-acceleration, and nQ{z) by a continuity 
equation as explained in the next paragraph. 



2.3 Pair production 

We compute at each altitude and frequency the pair pro- 
duction optical depth in the comoving frame. Absorption 
of 7-ray photons induces the formation of new pairs, that 
are supposed to be continuously reaccelerated as explained 
above. It results then in an increase of the particle flux 

$(2, ~ J ^i'y' ^' t)S{z)TbPbcd'y, which is conserved in the 

absence of pair creation. <l>(z) is computed through a con- 
tinuity equation : the importance of pair creation is mea- 
sured by the variation of this flux. The pair production op- 
tical depth is smaller than that expected in one-zone mod- 
els, where soft photons are produced cospatially with higher 
energy ones. Indeed a relativistic Maxwellian particle dis- 
tribution results in a much lower local soft photon density, 
compared to a power law spectrum : soft photons are pro- 
duced farther away in the jet and do not contribute to opac- 
ity. Consequently it allows signiflcantly lower bulk Lorentz 
factor compared to the ones obtained with simple one-zone 
models. In the case of PKS 2155-304 our results are com- 
patible with an asymptotic bulk Lorentz factor of about 15, 
much below the value of 50 infe rred by some authors (e.g. 
iBegelman. Fabian, fc Ree3l2008h . It turns out that pair pro- 
duction plays a fundamental role in explaining the flares 
in gamma-rays, because when the optical depth is close to 
one, a very small variation of the initial particle density can 
trigger a considerable enhancement in the pair production 
region. 



2.4 The jet spectrum 

Knowing 710(2, f), 70(2,4), r{z,t) and B{z,t), we compute 
the emissivity at each altitude in the jet by assuming a SSC 
process for the radiative mechanism. The total intensity of 
the jet is then determined by integrating the emissivity all 
along the jet. The emissivity is enhanced in the observer 
frame by the Doppler boosting. We take also into account 
the attenuation of the gamma ray signal by the cosmic dif- 
fuse infrare d background, chosen a s the "modified" Primack 
model P45 (|Aharonian et al.ll2006l '). 

Once injected at the base of the jet, the particles con- 
tribute first to the high energy part of the jet SED. As they 



propagate, their emissivity peaks progressively at lower en- 
ergy, producing the low energy part of the spectrum. Hence, 
the jet can be seen as a continuous succession of time depen- 
dent one-zone SSC models that propagate inside a station- 
ary geometry, each of them contributing with its particular 
spectrum to the whole observed SED. The spectral shape of 
the whole SED is not controlled by the local particle energy 
distribution (which is always locally a narrow pile-up), but 
rather by the z-dependencies of the jet radius, the magnetic 
field, and the acceleration rate. A constant injection rate 
would lead to a stationary emission pattern, which would 
be rather easy to fit. In reality however, the observed in- 
stantaneous spectra are a complicated convolution of the 
whole history of the jet, keeping the memory of the whole 
past (and unknown) injection pattern. We describe below a 
simple procedure to extract the physical parameters of the 
jet from observed spectra, despite the fact that they do not 
correspond to a simple steady-state of the jet. 

2.5 Time dependent simulations. 

Real observations result from a complex combination of the 
whole injection pattern at the basis of the jet. Due to the 
propagation of particles, the high energy part of the syn- 
chrotron and TeV components will be dominated by a sin- 
gle fiare, occuring at the basis, whereas the low energy part 
of the spectrum is a convolution over a large scale of the 
past jet history: it is thus rather a time-averaged spectral 
state mixing quiescent and fiaring states. The instantaneous 
spectrum is thus a mix of different ideal "steady-states". 

To constrain the geometrical parameters, we first build 
a virtual set of data called the fake fiaring spectrum, that 
would be observed if the jet were constantly fiaring. The 
high energy part of this spectrum is the really observed 
emission during a fiare. On the other side, its low energy 
part corresponds to the observed data corrected by an en- 
hancement factor to take into account that the actual jet 
is fiaring only a fraction / (called duty cycle) of the time. 
This enhancement factor is then since the real jet is 
filled only partially with the high particle density associated 
to a flare, whereas the fake spectrum corresponds to a fully 
filled jet. For the intermediate (X-rays) energy bands, the 
procedure is not so accurate because a small (but statisti- 
cally variable) number of flares can contribute to the flux. 
We chose to adjust approximately the spectra by some in- 
termediate factor in this band, but we have checked that our 
results are not very sensitive to this points. The / factor can 
be determined observationally from the fraction of time we 
see the object flaring within statistically random observation 



4 T. Boutelier et al. 



periods. Once the "fake flaring" spectrum is constructed, we 
use it to adjust the jet model geometrical parameters, and 
the flaring injection rate. We take also into account the re- 
quired variability timescale which constrains basically the 
jet size and the Doppler factor. The variability timescale is 
commonly evaluated in one-zone models by the light travel 
time of the zone divided by the Doppler factor R/Sc. How- 
ever, in the case of a stationary structure like a jet, the 
correct estimate is to take the typical length of the emis- 
sion zone , corrected for the time contraction due to the 
difference in light travel time to the observer along different 
parts of the jet, like for the superluminal motion. This gives 

Z Z 

tvar — — (1 — A cos 9) = which is of the same order of 

c cFftO 
the one-zone (transverse) variability timescale if we assume 
Ro/Zo ~ r^^. The model predicts also naturally an increas- 
ing variability with energy since more and more individual 
flares contribute to lower energy ranges. This is indeed ob- 
served in TeV blazars (jcicbcls ct al..,2007i ) . Once the flaring 
state parameters are found, we construct a quiescent state 
by keeping the same jet geometry, but reducing the par- 
ticle injection rate and/or acceleration rate to fit low flux 
observations. A general light curve can then be obtained by 
adjusting an injection pattern oscillating around the average 
flaring state, for a given particular flare. 
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Figure 2. Fit of PKS 2155-304 data. Filled dots: average 
archival data (see text). Empty triangles: average HESS data 
from the big flaring night. Empty diamonds: "fake flaring" state 
low energy points with a duty cycle f=0.1 (see text) Shaded 
area: enveloppe of archival X-ray data from BeppoSAX, SWIFT 
and XMM-Newton. Dot-dashed line: best fit of the " aver- 
age fake flaring" spectrum. Dashed line: flt of the quiescent 
(^average) spectrum. Solid line: example of an instantaneous 
simulated spectrum. 



3 APPLICATION TO PKS 2155-304 

We have tested our model to the big flar e event observed 
by H ESS in July 2006 in PKS 2155-304 (|Aharonian et all 
I2OO7I ). the strongest flaring event ever seen in the TeV range 
for a blazar. During this flare the average flux above 200 GeV 
of the source reached about ~ 7 crab. Moreover variabil- 
ity with timescales as small as ~200 seconds have been ob- 
served. The light curve is very well sampled (see Fig. |3]), 
as well as the spectrum, which makes this event particulary 
interesting to test models. 

We apply the method described previously to the 
SED of PKS 2155-304. We have first constructed an 
average spectrum of PKS 2155-304 in the low fre- 
quencies range (< 10^^ Hz), by compiling available 
archival data (from the HEASARC archive website 
I http : / /heasarc .gsfc .nasa.gov /docs/archive .html | ) . Based 
on the activity detected by HESS, we estimate that 
the duty-cycle is around 10%. An accurate value is not 
necessary since it would modify only the fake spectrum 
and not the real one. Then to construct the "fake flaring" 
spectrum of this source we combine the HESS data with the 
radio-to-optical ones corrected by a factor 10. To better 
constrain the space parameter we include also archival 
X-ray data from BeppoSAX, XMM-N ewton and SWIFT at 
different epoc hs dM assaro et al.' '20021), as well as archival 
EGRET data (|vS;trand ct al.. .1995 ). We choose not to 
re-normalize these data by any intermediate "duty-cycle" 
factor. Then they would correspond respectively to lower 
limits of the "flaring state" X-ray and GeV flux. 

An average "fake flaring" spectrum is shown in Fig. [2] 
in dot-dashed line. The corresponding best fit model pa- 
rameters are reported in Tab. [1] Interestingly we only 
needs a Doppler factor Vb ~ 15 which is significantly be- 
low the values of ~ 50 inferre d from one-zone model 
iBegelman. Fabian, fc Reesll200d ). For this simulation, the 



characteristic time scale Ro/{5i,c) is ~200 seconds in agree- 
ment with observations. 

Once we have obtained the best fit parameters of the 
"fiaring state" , we fit the TeV "quiescent state" by fixing the 
fit parameters associated with the jet structure and dynam- 
ics, allowing only the density of the particles filling the jet 
free to vary. We choose to fi t the averaged spectrum derived 
above 200 GeV by HESS in lAharonian et all (|2007l '). since it 
does not strongly differ from the quiescent state. The corre- 
sponding best fit parameters are also reported in Tab.[T]and 
the best fit model has been overplotted in Fig. [5] in dashed 
line. 



3.1 Time dependent simulation 

After having constrained the jet and plasma characteristics 
of PKS 2155-304 in both quiescent and fiaring states, we 
aim at reproducing the lightcurve observed by HESS during 
the big flare event by u sing a variable injection of particle. 
Based on the analysis of lAharonian et all (|2007h . where the 
light curve is decomposed in 5 successive assymetric bursts, 
we use, for the flare period, an injection function ^{z = 
0, t) that is the su m of five "generalized Gaussian" shape 
(|Norris et al.lll996l ) e.g.: 



<I>(0,t) = ^$^,cxp- 



(6) 



(jj the rise {t < t\^ax) a-nd decay {t > t^ai) time constant 
respectively, and is a measure of the sharpness of the 
burst. We have plotted in the middle of Fig. [3] the assumed 
injection function. Before the fiare, $(2; = 0,t)is assumed to 
be a crenel function that oscillates between quiescent and 
fiaring states in agreement with the source duty cycle. At 
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Figure 3. Upper panel: HESS light curve above 200 GcV su- 
perimposed with the model (solid line). Middle panel: time de- 
pendent particle injection function used in the simulation. Lower 
panel: predicted light curves in X-ray (dashed line, left y-scalc) 
and optical (dot-dashed line, right y-scale). The dotted lines mark 
the maximum of the different bursts of the injection function. 

the top of Fig. Owe have reported the HESS hght curve and 
the simulated one. The agreement is very good. We have also 
overplotted in Fig.[2]in solid line an instantaneous spectrum 
extracted from the time-dependent simulation during the 
flare period. It agrees nicely with the broad-band (from radio 
to TeV) spectrum observed during this burst. 



4 DISCUSSION 

Our time-dependent inhomogeneous jet model succeeds in 
reproducing simultaneously the broad band (from radio to 
TeV) spectrum of PKS 2155-304 as well as the TeV Ught 
curve during the big flare event of July 2006. The key idea 
of the method is to decompose the blazar spectrum in "qui- 
escent" , low luminosity states, and "flaring" , high luminosity 
states. The high energy part of the spectrum, coming from 
small-scale inner regions, is assumed to be, at any time, in 
one of these pure states. On the other hand the low energy 
part is a convolution over a large scale of the past history 
of the jet : it is thus rather a time-averaged spectral state 
mixing quiescent and flaring states in proportions given by 
the source "duty-cycle". Moreover we do not require two 
differen t popu lations of emitting part icles like in other m od- 
els (e.g. iKatarz viiski ct al. 2003, Cl iiaberge et al.ll2000l ) but 
simply a continuous (although variable) injection of parti- 
cles at the base of the jet that propagate along the jet struc- 
ture. Pair production plays an important role to amplify the 
initial variation, as can be seen with the variation of the par- 
ticle flux along the jet during the flaring state (see Tab. 1) 
: the pair current is amplified by a factor 30 at the end of 
the jet, when the initial current varies only by a factor 2. 

The model can also predict light curves at different 
wavelengths. As an example, the X-ray (2-10 keV) and op- 



tical (V band) light curves expected during the TeV flare 
have been plotted at the bottom of Fig. [3] The X-ray lu- 
minosity exhibits almost simultaneous variations but with a 
lower amplitude (about 5 times smaller). On the other hand, 
the optical light curve shows a very different behavior, in- 
creasing all along the flare. This is due to to the large size 
of the optical emitting region that plays the role of a low 
pass filter. Consequently, the optical luminosity integrates 
the recent past history of the jet. These results are compat- 
ible with simultaneous multiwavelength observations made 
during the "Chandra night" (oral communication). 

This model puts also some constraints on the central 
black hole mass. Indeed, the radius at the base of the jet 
is directly linked to the bulk Doppler factor. From our 
best fit parameters, we find Mth = 8 x lO^Vg/Ro, where 
Tg = GMiyh/c^ is the gravitational radius. This implies a 
black hole mass smaller than 10^ Mp, usually in voked in the 
case of PKS 2155-304 (|Kotilainen et al.lll998l ). We believe 
however that it is s till acceptable given the mass measure- 
ment uncertainties (jBettoni et al.ll2003l ). 

In the present work, we do not specify the origin of the 
variability. Obviously, other plasma parameters can vary in 
addition to the injection density A'^o. Most likely, variabil- 
ity can be triggered by a change in the acceleration rate 
described by Qq (Eq. O. In a plausible scenario, long term 
(year scale) variability implying the succession of quiescent 
and active states could be attributed to variations in the ac- 
cretion rate, wheras the short (minute-scale) fiares would be 
attributed to the instability to pair creation that develops 
only when the initial particle density is close to a critical 
threshold. 
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